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Objectives: Bacteroides fragilis, a Grann-negatlve anaerobic bacterium, is alternately a gut connmensal or virulent 
pathogen and is an important reservoir for horizontal gene transfer (HGT) of bacterial resistance and virulence genes in 
the human gastrointestinal tract. We identified a unique conjugative transposon (CTn) in a multidrug resistant clinical iso- 
late of 6. fragilis (BF-HMW615); we named this element CTnHyb because it included a hybrid mosaic of foreign elements. 
This study reports the characterization of CTnHyb and discusses the potential impact on horizontal spread of resistance 
genes. 

Results: CTnHyb contains several efflux pump genes and several genes that confer or may confer antibiotic resis- 
tance to tetracycline, kanamycin, metronidazole and spectinomycin (truncated gene). CTnHyb also contains a mosaic of 
mobile elements from Gram-positive organisms. CTnHyb is easily transferred from BF-HI\/1W615 (the original isolate) to 
BF638R (lab strain) and integrated into the BF638R chromosome. The "foreign" (from Gram-positive bacteria) nucleotide 
sequences within CTnHyb were > 99% preserved indicating that the gene acquisition from the Gram-positive bacteria 
was very recent. 

Conclusion: CTnHyb is a novel CTn residing in a multidrug resistant strain of 6. fragilis. The global nature and wide 
phylogenetic reach of HGT means that any gene in any bacterium can potentially be mobilized. Understanding the 
mechanisms that drive the formation and transfer of these elements and, potentially, ways to limitthe transfer are neces- 
sary to prevent a devastating spread of resistance elements. 



Introduction 

The fluidity of the human gut microbiome has been recog- 
nized for decades but the recent data explosion from the analysis 
and sequencing of the gut microbiota is clarifying the vast extent 
of this gene transfer. Not only are genes regularly transferred 
among permanent residents of the gut, but organisms such as 
Staphylococcus and Streptococcus that do not colonize the intes- 
tine can participate in this genetic "swap meet"' (when they pass 
through the gut). 

Of the resident gut population, -99% of the microbiota 
belongs to two divisions (superkingdoms) of Bacteria-the 
Bacteroidetes (48%) (including Bacteroides species) and the 
Firmicutes (51%).' Bacteroides species, the dominant bacterial 
genus in the human gut are known to harbor many conjugative 
and mobilizable elements.^ 

Conjugative transposons, also known as integrative conjuga- 
tive elements (ICEs), are a subset of mobile elements that also 



include plasmids and transposons.^ Like transposons, CTns can 
integrate into diverse sites in the host chromosome. The CTns 
do not exclude each other as do plasmids, so a strain can accu- 
mulate more than one CTn. Furthermore, there is some evi- 
dence that the presence of more than one copy of the CTn in 
the strain results in a stimulation of transposition (transactiva- 
tion).'' Theoretically, this implies that as CTns with antibiotic 
resistance genes accumulate in the environment, the transfer of 
these genes to other bacteria will also increase and may result in 
upward spiraling of antibiotic resistance.^ 

We recently investigated a multidrug resistant clinical isolate 
of Bacteroides fragilis (BF-HMW615) isolated from a pediatric 
appendiceal specimen.'' BF-HMW615 is resistant to multiple 
antibiotics, including metronidazole. Our analysis identified 
a nimj gene which conferred increased MICs to metronidazole 
when introduced into a susceptible strain^ and nimJ is carried 
on a novel conjugative transposon. We are now reporting the 
identification of this novel conjugative transposon, CTnHyb (for 
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Figure 1. Genomic sequence comparison using the RAST Server (http://rast.nmpdr.org/). HiVIW 615 is compared with HMW 610, HIVIW 616, BF 9343, 
BF638R, and BF YC[H46, respectively from perimeter of circle inward; HMW 615 sequence is inferred (not shown). The sequence of HIVIW 615 has been 
taken from the published supercontigs and rearranged according to the Mauve prediction. Hits on the comparison organisms are displayed graphically. 
Percent protein sequence identity is indicated (legend). 



"hybrid"), which contains genes from Gram-positive bacteria. 
CTnHyb is transferable to B. fragilis 638R and thus is confirmed 
as a mobile element. The CTnHyb has (besides from Bacteroides 
spp) exact nucleotide homologs from at least three phyloge- 
netically distinct Gram-positive organisms. The extent of the 
"hybrid" nature of this CTn, to our knowledge, has not been 
reported before in Bacteroides. 

Results 

Identification of CTnHyb in BF-HMW615 by comparative 
genome analysis 

Comparative RAST-based genomic analysis'*'' indicated a 
few continuous regions of BF-HMW615 chromosome (> 50,000 
bp) that contained genes with no homologs in B. fragilis ATCC 
9343, BF638R, BF Y46H or the multidrug resistant clinical iso- 
lates BF-HMW610 and BF-HMW616. The blue lines indicate 
that BF-HMW615 has homologs and red lines indicate absence of 
homologs (Fig. 1). One of these segments (indicated by red arrows. 
Fig. 1) included the tetCl gene as well as multiple Bacteroides tra 
genes (implicated in conjugative transposition). Thus, we con- 
sidered this a potential CTn. Initial BLAST analysis of this non- 
homologous segment in BF-HMW 615 indicated that some of the 
genes have been horizontally transferred from other species. 

Transfer of CTnHyb to B. fragilis 638R 

We mobilized CTnHyb from BF-HMW615 to BF638R by 
mating. The new 638R-CTnHyb mutant was selected using 



tetracycline and rifampicin (BF638R is rifampicin resistant, 
tetracycline sensitive; BF-HMW615 is rifampicin sensitive, tet- 
racycline resistant). To rule out the small chance of the tetR/ 
rifR transconjugant of BF-HMW615 origin, we confirmed the 
BF638R origin by PGR amplification and partial sequencing of 
two BF638R genes (BF638R_2089 and BF638R_4382) which 
are not present in BF-HMW615. There was no evident tetracy- 
cline-mediated increase in frequency of CTnHyb. 

Determination of the CTnHyb ends and the point of inser- 
tion into the BF638R chromosome 

We selected six colonies (HMW 874, HMW 875, HMW 
876, HMW 877, HMW 878, and HMW 879) from indepen- 
dent mating experiments (of BF638R and BF-HMW615) to 
determine the insertion point of CTnHyb. The insertion points 
were determined by the SRP technique as described in Materials 
and Methods. The isolates HMW 874, HMW 875, HMW 877, 
and HMW 878 had identical ends and insertion points, whereas 
HMW 876 and HMW 879 had slightly altered ends and inser- 
tion point. 

The genome of BF-HMW615 has been sequenced but the 
published sequence is in supercontigs: the sizes range from 1.1 
(largest supercontig) to 1.14 (smallest supercontig). Our analysis 
indicated that the transferred CTnHyb fragment was partially 
contained in two supercontigs annotated as 1.1 and 1.3. Our 
PGR data connecting the supercontigs 1.1 and 1.3 was supported 
by MAUVE alignment of BF-HMW615 supercontigs with 
BF638R (Fig. 2). Also, an unannotated transposase was identi- 
fied, through PGR, in the junction sequence of 1.1 and 1.3. 
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Figure 2. MAUVE-based alignment of BF 638R and BF-HMW615 and schematic indication of CTnHyb insertion into BF638R. The MAUVE Contig IVlover 
module was used to order (and orient) the BF-HMW615 supercontigs relative to the BF 638 reference genome. Colored block is a region of the genome 
that is homologous to a similarly colored part of another genome without extensive genomic rearrangement. Regions outside blocks lack detectable 
homology among the input genomes. Only a small part of the alignment is shown to visualize the insertion; n/mJ and fetQare indicated for orientation. 



The predicted crossovers that would result in the sequences in 
638R/CTnHyb isolates is depicted in Figures 3A and B, respec- 
tively. The left end of CTnHyb is dnaK2 (HMPREF1204_0027 
in reverse complement supercontig 1.1 of BF-HMW615) 
and the right end is between HMPREF1204_03040 and 
HMPREF1204_03041 (in contig 1.3 in BF-HMW615). The 
insertion is in dnaK2 (BF638R_1256, between 1507197 and 
1507222 bp). 

Determination of size of CTnHyb 

The two transconjugants, HMW 874 and HMW 875, were 
used to determine the length of CtnHyb. Surprisingly, two dif- 
ferent size inserts had been incorporated into the BF638R chro- 
mosome. HMW 875 had a larger insert (thus named, CTnHybL) 
of 131,471 bp and HMW 874 had a smaller insert (thus named, 
CTnHybS) of 98,099 bp. Sequence analysis indicated that both 
inserts had the same ends but there were two internal segments 
missing in the CTnHybS (98099 bp), leading to the transconju- 
gant with the smaller insertion as depicted in Figure 3C. Primers 
were designed across distances of 13327 bp (if the first deletion 
had not occurred) and 21512 bp (if the second deletion had not 
occurred), respectively, to detect whether the deletions were pres- 
ent in various isolates. PCR and subsequent sequencing across 
these regions indicated that the deletions were present not only 
in HMW874, but were also present in chromosomal DNA of 
BF-HMW 615 and in HMW875. The amplicons spanning the 
deletions (861 and 579 bp, respectively) were sequenced to con- 
firm that they represented the deletion region. Thus, it appears 
that these deletions can occur in BF-HMW 615 as well as in 
BF638R/ CTnHybL indicating the volatile nature of CTnHyb. 
Determination of the circular form of CtnHyb 
The primers (2783 and 2784) were used in a PCR reaction 
on BF-HMW6I5 and the BF638R/CTnHyb transconjugants 
to determine if we could detect a circular intermediate, which 
is a critical intermediary step in CTn transfer. The resultant 



amplicons indicated that there was indeed a junction between 
the two ends of CTnHyb indicating the presence of circu- 
lar intermediate in BF-HMW 615. Interestingly, the circular 
intermediates were also detected in 4 independent isolates of 
BF638R/CTnHyb which indicated that CTnHyb can continue 
to make conjugative circle intermediates even after transfer to 
BF638R. 

To visualize this circular intermediate, we constructed a 
GenBank (gb) file for CTnHyb from: 1) the sequence included 
within supercontig 1.1, 2) an unannotated transpose downstream 
of the predicted HMPREF1204_0001, 3) a short sequence 
obtained by PCR of the supercontig 1.1 and 1.3 junctions, and 
4) the sequence contained within supercontig 1.3. The 4 gb files 
were concatenated using the SeqNinja program (DNASTAR, 
Inc., Madison, WI). The concatenated file was annotated and 
viewed using SeqBuilder (DNASTAR, Inc., Madison, WI) 
(Fig. 4). The nimJ and tetQ genes are indicated, as are the "for- 
eign element" (from Gram-positive bacteria) extending from 
HMPREF1204_2965 to HMPREF1204_2980 and the pre- 
dicted crossover point. The inner most circle is the GC% (green 
is higher than average, purple is lower than average) . 

CTnHyb contains 144 genes 

BLAST analysis of the CTnHyb indicated that it contained 
CTn specific tra genes {traE, traG, traj, traK, traM and tra^, 
excisionases, transposases and other DNA-associated proteins, a 
tetracycline resistance gene (tetQ), 3 putative pump system genes 
coding for efflux pumps (MefA, ABC, and RND type transport- 
ers), genes coding for hemagglutinin and thioredoxin (both may 
be important in virulence), and genes coding for metronidazole, 
kanamycin and tetracycline resistance. The genes contained 
within CTnHyb are shown in Table 52. RNA-Seq results for all 
of the BF-HMW 615 genes will be published as part of a larger 
study comparing the total transcriptome of BF clinical isolates 
(Husain F, Veeranagouda Y, Wexler HM, unpublished data) but 
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Figure 3. (A) Predicted events leading to integration of CTnHyb into the BF638R chromosome to result in BF638R/CTnHyb. The circular form of CTnHyb 
recombines with BF638R at the underlined sequence (GAAAGTGAA). (B) Alternate integration of CTnHyb into the BF638R chromosome. Predicted 
events leading to integration of CTnHyb into the BF638R chromosome to result in BF638R/CTnHyb found in HIVIW 876 and HMW 879 . The circular form 
of CTnHyb recombines with BF638R at the underlined sequence (TTTGA). (C) Predicted model of the deletions in CTnHybL leading to CTnHybS. The 
locus tag labels serve as approximate reference points. The bases are referred to as "bases" or "b". The regions of homology that are predicted to recom- 
bine are color coded. The nimJ gene is among the genes deleted in the first deletion and metronidazole has a lower MIC for BF638R/CTnHybS than for 
BF638R/CTnHybL, as expected (Table 2). 



the RNA-Seq counts for the CTnHyb genes are presented here, 
as well. The distribution of genes in CTnHyb according to COG 
(Cluster of Orthologous Gene class) is shown in Figure 5. 

CtnHyb includes a "foreign segment" containing genes 
homologous to a variety of Gram-positive bacteria (Fig. 4; 
Fig. 6A). HMPREF 1204_02969 (coding for amino- 
glycoside 3 '-phosphotransferase conferring kanamycin 



resistance) is 100% identical to Staphylococcus epidermidis 
RP62A aphA (SEAOOlO). HMPREF 1204_02965 encodes 
a MefA type efflux pump and is 100% homologous to mefA 
genes in many Streptococcus species. Downstream of these 
genes is a 6790 bp segment (HMPREF1204_2967 through 
HMPREF1204_02977) that is homologous to a nucleotide 
segment in Eubacterium ventriosum 27560 (a Gram-positive 
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Figure 4. CTnHyb: Circularized form and predicted nested ICEs. Circular intermediates were detected in HMW615, and in 4 isolates of BF638R/CTnHybL 
using CTnHyb/BF 638R Junction primers at either end of CTnHyb (Primers 2783 and 2784). The continuous sequence was generated using 4 GenBank 
files consisting of 1) the known Broad sequence within supercontig 1.1, an unannotated transpose downstream of the predicted HIVIPREF1204_0001, 
a short nucleotide segment and the sequence within supercontig 1 .3 (we obtained the 2nd and 3rd sequence from manual PCR and sequencing). The 
GenBank files were concatenated using the SeqNinja program (DNAStar) and the generated sequence was annotated and visualized with the SeqBuilder 
program (DNA Star). The inner most circle is the GC% (green is higher than average BP GC % and purple is lower than average BP GC %). The predicted 
attachment site, the nimJ, mefA and tetQ genes and the predicted nested ICEs are indicated on the figure. 1) CTnHyb "Popout 1": (missing in CTnHybS) 
12466 bases long, extending from HIV1PREP1204„00014 to HMPREP1204_00001 of contig 1.1 and the unannotated transposase and HMPREP1204_02924 
to HIV1PREP1204_02926 of contig 1.3); 2) NimJ cassette, 9955 bases, present in 3 copies in BP HIVIW 615; 3) Cassette (minus nimJ) is found in other BP 
strains; 3) nimJ gene; 4) CTnHyb "Popout 2 (missing in CTnHybS): 20,933 bases long; extending from HMPREP1204_02928 to HIVIPREP1204_02946; 5) 
The "foreign segment" extends from HMPREF1204_2965 to HI\/1PREP1204_2980; 6) mefA (HIVIPREP 1204_02965, 1232 bp) is a 95% match to mefA genes 
on ICEs from Staphylococcus, Streptococcus, and Enterococcus^-^, and a close mefA homolog was also found in CTnGerm^'; 7) Eubacterium ventriosum cas- 
sette (7178 bp not counting the kanamycin cassette within); 8) kanamycin cassette (1517 bp); HMPREP1204_2969-HMPREP1204_2971 is 100% match to 
nucleotide regions on ICEs from Streptococcus. However, the transposase (HIVIPREP 1204_02971) in the "kanamycin cassette" is not homologous to the 
Staphylococcus or Enterococcus transposases, but is homologous to a transposase from Blautia hansenii — a novel genus of Gram-positive, anaerobic, 
non-sporulating coccobacillus-shaped bacteria that includes several former coccobacillary shaped species of Clostridia and Ruminococcus.-"' HMPREP 
1204_02971 is truncated at a 10 bp palindrome (ACTTCCGCCG, bp 195-204 and 218-227 within HMPREP 1204_2970) which is characteristic of trans- 
posase insertions. Bacterial repetitive extragenic palindromic sequences are known DNA targets for Insertion Sequence elements. This palindrome, 
however, is contained within the coding region of HMPREF1204_02971. 



gut anaerobe) from EUBVEN_02875 thru EUBVEN_02862 
except for a short missing stretch (Fig. 6B). The "missing 
piece" is replaced by a 3617 bp sequence (HMP1204_02969, 
HMP1204_02970 and HMP1204_02971, GenBank Accession 
AGXR01000023.1 46489-48006) that is highly con- 
served in both genome sequences and ICEs in Staphylococcus, 
Streptococcus and Enterococcus. In S. pneumonia, for example, 
it is present on mobile elements carrying multidrug resistance 
determinants''' and in Enterococcus faecalis RE25 it is present 
on a 50-kb conjugative multidrug resistance plasmid (pRE25)'^ 
(Fig. 6B). We don't know whether the insertion of the 3617 
bp segment (containing HMPREF1204_02969, 02970 and 
02971) happened before or after the Eubacterium segment was 
transferred into BF-HMW 615. Also, whether the acquisition of 
the E. ventriosum genes is a result of a BE CTnHyb-like element 



having moved into Eubacterium and then back into BE is not 
clear at this point. 

CTnHybL contains a novel metronidazole resistance gene 

CTnHyb contains a cassette with a recently reported metro- 
nidazole resistance ("«/w/") gene, found in two metronidazole 
resistant clinical BF isolates (BF-HMW 615 and BF-HMW 
616).^ Downstream of nim] is an unannotated transposase of 
the IS4 class (in opposite orientation) which spans the junction 
of supercontigs 1.1 and 1.3. It is in reverse orientation to nim] 
and contains the consensus promoter sequence TAnnTTTG that 
is found on insertion sequences containing the cfiA (imipen- 
emase) gene."" nim] (HMPREF1204_00002) is present within 
CTnHybL but is not present in CTnHybS (see Fig. 3C). As 
expected, BF638R/CTnHybL had increased MICs for metro- 
nidazole (1.0-1.5 JJtg/ml) while BF638R/CTnHybS had MICs 
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Figure 5. Distribution of genes in CTnHyb according to COG (Cluster of Orthologous Gene class). COG ID numbers and classification were taken from the 
Integrated Microbial Genomes and Metagenomes (IMG) at the Department of Energy Joint Genome Institute, (http://img.jgi.doe.g0v/cgi-bin/w/main. 
cgi?section=FindGenes&page=geneSearch). Eighty-eight genes in CTnHyb do not belong to any COG. Similarly, R AST analysis of the entire BF-HMW615 
genome sequence indicates that 69% of the putative genes are not assigned to any subsystem (a subsystem is a method of categorization that can be 
thought of as roughly equivalent to a COG). Of the remaining 57 genes in CTnHyb(L), almost one third of the genes (n = 18) are in COG L (replication, 
recombination and repair) which includes such genes as DNA primases, excionases, integrases and transferases. Another 11 genes are annotated as 
"viral proteins" and are not assigned to a COG class, even though several are also functionally annotated as excisionase or integrase proteins. 



Figure 6 {See opposite page). (A) Schematic representation of the "foreign" segment in CTnHyb. White arrows correspond to conserved genes in 
Bacteroides sp strains. See the figure for explanation of other arrows. APH(3'): aminoglycoside 3'-phosphotransferase; ABC: ABC transporter ; Mob: 
mobilization protein; Spe R: spectinomycin adenylyltransferase; KNTase: nucleotidyltransferase; HMTase: methlytransferase; HDc: Metal dependent 
phosphohydrolases with conserved "HD" motif; SIS: Sugar Isomerase; ZnCP: Zinc peptidase; BmgA and BmgB: mobilization proteins. (B) ACT alignment 
of area of homology between Eubacterium ventriosum and CTnHyb (BF-HMW615). A segment (7671 bp) of the E. ventriosum genome is nearly completely 
conserved except for a three gene insert that replaces part of EUBVEN_02872, EUBVEN_02873 and part of EUBVEN_02874. In the conserved portions, 
only 11 nucleotides differed. 



Table 1. Minimal inhibitory concentrations of strains and constructs 





MICs (mg/l)» 




Tetracycine 


Erythromycin 


Metronidazole 


BF 25285 


0.60 ± 0.22 


3.5 ±0.71 


0.41 ± 0.04 


BF 638R 


0.16 ±0.00 


1 ± 0.00 


0.44 ± 0.08 


BF HMW615 


192.00 ±0.00 


256 ±0.00 


12 ±5.66 


BF 638R/CTN HybL 


60.00 ± 5.66 


13 ±4.24 


1.125±0.18 


BF 638R/CTN HybS 


64.00 ± 0.00 


6 ± 0.00 


0.41 ± 0.04 



'MIC is the average of two technical repeats and two biological repeats. 
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in the same range as BF638R (0.38-0.5 |xg/ml). Additionally, 
MICs of erythromycin for BF638R/CTnHybL and BF638R/ 
CTnHybS were increased from those of BF638R (Table 1). 

CTnHyb contains multiple transposase genes of the IS4 
category 

In addition to the transposase gene adjacent to nimj, 
CTnHyb contains two additional transposase genes of the IS4 
class that are divergently transcribed from their adjacent genes. 
HMPREF1204_2983 is adjacent to HMPREF1204_2984 
(tetQ) and contains the consensus promoter sequence, and 
HMPREF1204_2964 is another IS4 family transposase adjacent 
to HMPREF1204_2965 (mefA). 

Expression of HMPREF1204_02969 (conferring kan") and 
HMPREF1204_02965 (mefA) in E. coli 

HMPREF1204_2969 (encoding aminoglycoside 

3 '-phosphotransferase for kanamycin resistance) and 
HMPREF1204_02965 {mefA) were cloned into pSportI 
and introduced into E. coli: HMPREF1204_2969 into E. 
coli AGIOO to determine changes in kanamycin MICs and 
HMPREF1204_02965 into both E. coli AGIOO and E. coli 
Kam43 (a pump deficient mutant, lacking AcrAB, AcrEF 
and TolC'^) to determine erythromycin MIC changes due 
to mefA expression. E. co/z7pSportI::HMPREF1204_2969 
grew in the presence of 40 [xg/ml kanamycin, while E. coli! 
pSportI (without the kanamycin gene) did not grow, indi- 
cating that HMPREF1204_2969 was fully functional in E. 
coli. HMPREF1204_02965 {mefA) did not confer erythro- 
mycin resistance in either AGIOO or KAM43. The mating of 
BF-HMW615 and E. co/z AGIOO did result in DNA transfer that 
yielded E. coli that were kanamycin resistant, but the transcon- 
jugants were not stable and did not yield further generations on 
purification. 

Discussion 

The importance of the gut flora in human health and dis- 
ease is at the forefront of scientific and public awareness and 
major efforts are underway to characterize and sequence the 
gut microbiota as part of the massive Fluman Microbiome 
Project (HMP) at the Broad Institute/NIFI. The gut bacte- 
ria, in terms of total cell and gene numbers in the human"*, 
"represents a virtual inner organ"". BE is a human gut com- 
mensal that only accounts for 2% of the total Bacteroides but 
it is the agent of > 70% of Bacteroides infections.' As a com- 
mensal, it hydrolyzes complex polysaccharides and produces 
volatile fatty acids used by the host as source of energy' and 
is important in immune development. However, it is very 
virulent when it escapes the gut and has been associated with 
nearly all types of infections. Additionally, it may be a res- 
ervoir of resistance genes that can get passed, by horizontal 
gene transfer (HGT), to other organisms resident in or passing 
through the gut. 

Horizontal gene transfer among gut microbiota is particu- 
larly intense; gut microbes, therefore, may be a major reservoir 
for antibiotic resistance genes. Indeed, the taxonomically 



different representatives of gut microbiota may share the pool of 
closely related antimicrobial resistance genes. HGT is also a cru- 
cial event in the development of virulence traits.^'' 

One of the major elements responsible for HGT are conjuga- 
tive transposons (CTns). CTns are similar to transposons (that 
integrate into the host semi-randomly) but also carry the genes 
necessary for conjugal transfer to other cells. The first complete 
sequence of the transfer region of a Bacteroides conjugative trans- 
poson was described in 2001.^' The many known conjugative 
transposons and other mobile elements that are present in strains 
of Bacteroides were reviewed recently^ and include CTnDOT 
(the most widely studied Bacteroides CTn),""^ CTnGerm," 
CTn34l,^^ CTnBST-« and CTnl2256.-' CTnGerm and CTnBst 
were reported to carry a variety of genes with high similarity to 
genes from aerobic bacteria"'^* and CTnl2256 was described as 
a chimeric transposon composed of two independently active 
mobile elements.^' Sequence comparisons between CTnHyb 
and other Bacteroides CTns, including CTnDOT, CTnGerm or 
CTnBST did not reveal any significant homology. 

CTnHyb is unique and noteworthy among the many conju- 
gative transposons and other mobile elements present in strains 
of Bacteroides. First, several genes coding for antibiotic resis- 
tance proteins are contained within CTnHyb. The tetCl and 
nimJ genes conferred tetracycline and metronidazole resistance, 
respectively, in BF638R and HMPREF1204_2969 conferred 
kanamycin resistance in E.coli. The encoded efflux pumps 
including the RND and mefA pumps may also contribute to 
drug resistance. Also, there is a truncated gene for spectino- 
mycin resistance. Second, CTnHyb has a long sequence stretch 
of homologous sequences from a Gram-positive species (i.e., 
Eubacterium) that includes a mosaic of resistance genes from 
other aerobic Gram-positive species. The degree of nucleotide 
conservation (6779/6790 nucleic acids) within the area homolo- 
gous to Eubacterium ventriosum indicates a very recent transfer 
to BE, otherwise some degree of adaptation to codon usage bias 
patterns of BE would be expected.'" The manner in which the 
various "foreign" elements are arranged within CTnHyb sug- 
gests an "Ellis Island" effect where incoming CTns are prefer- 
entially drawn to genomic regions that are composed of other 
mobilizable elements; ICEs are nested within each other (Fig. 4) 
(note: Ellis Island was the gateway for millions of immigrants 
to the United States from 1892 to 1954). This type of modu- 
lar transfer within conjugative elements has been described in 
other bacteria and in one CTn in Bacteroides and may be an 
important mechanism for the accumulation of resistance and 
virulence genes in the gut and the subsequent development of 
resistance or pathogenicity islands. '' 

Target site selection is an important characteristic for each 
transposon and determines dissemination and stability.'^ We 
detected a circular intermediate of CTnHyb in BE HMW 615 
and in HMW615/BF638R transconjugants and predicted two 
potential crossover points that would yield sequences consistent 
with the sequences found in the resultant BF638R/CTnHyb 
transconjugants. In CTnDOT, the most extensively studied 
Bacteroides CTn, a tyrosine recombinase called IntDOT cata- 
lyzes integration into, and excision out of the bacterial host 
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chromosome. The core (GTANNTTT), are 
inverted repeat sequences that flank target 
sites in the chromosome and in CTnDOT, 
where strand exchange takes place catalyzed 
by IntDOT. The target sites, attB, on the host 
chromosome consist of a pair of inverted repeat 
core sites (B and B'). The complementary sites 
on CTnDOT, attDOT sites, have the core 
sites D and D'. CTnBST, another Bacteroides 
CTn, appears to integrate more site specifically 
than CTnDOT, with a 6-amino-acid signature 
that is associated with the catalytic regions of 
members of the tyrosine recombinase fam- 
ily.^- Although the core site sequence is present 
('5037 times in BF638R), it is not at the cross- 
over position of CTnHyb. 

Transfer of several of the widely studied 
Bacteroides CTns, including CTnDOT, is 
mediated by low levels of tetracycline and it is 
believed that the wide use of antibiotics thera- 
peutically and in animal feed (with subsequent 
contamination of both meat and manure- 
fertilized crops) influence the introduction of 
these mobile-element bearing-organisms into 
the human gut. It is possible that elimina- 
tion of the inducers (i.e., antibiotic) might be 
achieved by radically reducing their use. In the 
case of CTnHyb, however, transfer (from the 
multidrug resistant clinical isolate to the sus- 
ceptible lab strain) could be easily achieved without any need 
for induction by tetracycline so change in tetracycline use may 
not affect transfer frequency. The global nature and wide phylo- 
genetic pool of the horizontal transfer described in recent years 
means that any gene in any bacterium can potentially be mobi- 
lized and resistance phenotypes can be established in a diverse 
range of organisms worldwide. Understanding the mechanisms 
that drive this transfer and ways to limit the transfer are neces- 
sary to quell the spread of resistance elements. ^""''^'^^ 



Materials and Methods 

Strains and culture conditions 

Strains used in this study are listed in Table 2. All strains 
were grown as described'"' using Brain Heart Infusion media 
supplemented with 15 iJLg /ml hemin (BHIS) for Bacteroides iso- 
lates (Anaerobe Systems, Morgan Hill, CA) and Luria Bertani 
(LB) agar or broth (Sigma) for Escherichia coli. The multi- 
drug resistant clinical isolates BF-HMW610, BF-HMW615 
and BF-HMW616 have been previously described.*'^''''* E. coli 
AGIOO was used as the host to test for the kanamycin {kan) and 
spectinomycin resistance phenotypes. E. coli Kam43 (a pump 
deficient mutant, lacking AcrAB, AcrEF and TolC'^) was a kind 
gift from Dr Tomofuso Tsuchiya (Okayama University, Japan). 
Ampicillin (50 jJLg/ml), erythromycin (10 |JLg/ml), and kanamy- 
cin (40 (JLg/ml) were used for selection as indicated. 



Table 2. Strains and plasmids used in this study 



Strains 


Description 


Phenotypes 


Source or reference 


B. fragilis 








BF638R 


Wild type parental strain 


Tet^ Rif 


37 


HIV1W615 


Drug resistant clinical isolate 


Tet«, Rif 


6 


HIVIW 874 


638R/CTnHybS 


Tet«, Rif 


This study 


HIVIW 875 


638R/CTnHybL 


Tet", Met", Rif" 


This study 


HIVIW 876 


638R/CTnHyb 599 


Tet", Met", Rif" 


This study 


HMW 877 


638R/CTnHyb 600 


Tet", Met", Rif" 


This study 


HMW 878 


638R/CTnHyb 604 


Tet", Met", Rif 


This study 


HMW 879 


638R/CTnHyb 605 


Tet", Met", Rif 


This study 


E. coli 








AGIOO 


Cloning strain 


Amp^ Kan^ 


38 


DH5a 


Cloning strain 


Amp^ Kan^ 


Invitrogen 


Kam43 


TG 1 hacrAB, AydhE, AtolC, 
efflux pump deficient 


Amp^ Kan^ 


39 




AGIOO/CTnHybL 


Kan" (unstable) 


This study 




DH5a/pSP0RT1 -.-.kan 


Amp", Kan" 


This study 




DH5a/pSP0RT1 "mefA 


Amp" 


This study 




AG100/pSPORT::meM 


Amp" 


This study 




Kam43/pSPORT::meM 


Amp" 


This study 


Plasmids 


pSPORTI 




Invitrogen 




pSport\::kan 




This study 




pSportl::meM 




This study 



Molecular methods 

DNA extraction, restriction digestions, gel electrophoresis 
and analysis were done as previously described.''" The size and 
sequence of the transferred CTn was determined by PCR and 
sequencing of the CTn at regular intervals. Based on BF-HMW 
615 sequence, primers were designed to yield 150 to 200 bases 
products targeting DNA approximately every -20 KB on either 
side of the tetQ gene (HMPREF1204_02983). Using these 
primers, the PCR amplification was done with genomic DNA 
from BF-HMW 615, BF638R, and the selected BF-HMW 615/ 
BF638R transconjugants as templates. The exact boundar- 
ies of CTnHyb insertion into BF638R, the insertion points in 
BF638R, the segments deleted in CTnHybS, and the sequence 
of the junction of contig 1.1, the unannotated transposase, and 
contig 1.3 were determined by semi-random priming (SRP)- 
PCR.''^ The primers used for sequencing the BF638R/CTnHyb 
junctions and the gaps found within the CTnHyb in one of the 
isolates are listed in Table SI. 

Genome sequencing 

BF-HMW615 (along with two other multidrug resistant iso- 
lates BF-HMW610 and BF-HMW616) were submitted to the 
Broad Institute and sequenced as part of the Human Microbiome 
Project, Bacteroides group Sequencing Project, Broad Institute 
of Harvard and MIT (http://www.broadinstitute.org/). The 
Broad sequencing project utilized 454 Whole Genome Shotgun 
methodology and Newbler (454 Life Sciences) assembly. This 
sequencing project was supported by the National Institute of 
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Allergy and Infectious Disease/National Institutes of Health- 
funded Genome Sequencing Center for Infectious Diseases at the 
Broad Institute. BF-HMW610, BF-HMW615 and BF-HMW616 
have been given the Broad designations HMPREF1203, 
HMPREF1204, and HMPREF1205, respectively. For the sake 
of consistency, the BF-HMW 615 genes are referred to by their 
designation 1-IMPREF1204_, etc.; these are the designations 
used in GenBank FASTA files of the genome sequences and asso- 
ciated annotations were downloaded from the Broad Institute. 
Genomic analysis 

The RAST (Rapid Annotation using Subsystem Technology) 
Annotation Server** was used for comparative genome analysis. 
All sequences submitted to RAST were downloaded from either 
the Broad Institute (for the clinical isolates) or NCBI (for the ref- 
erence strains) and submitted to the RAST server.^' Genomes 
were compared by the sequence comparison feature of the 
SEED server.' Comparisons of known CTns and the BF-HMW 
615 genome sequence were analyzed as described earlier.^ The 
Double ACT server (http://www.hpa-bioinfotools.org.uk) was 
used to generate the comparison file with the BLASTN func- 
tion and a cutoff of 1000-2000. Results were viewed using the 
ACT viewer.'"'' MAUVE alignment software, which is specifically 
designed for multiple genome alignment in the presence of large- 
scale evolutionary events such as rearrangement and inversion,''^ 
was used to help predict the arrangement and orientation of 
supercontigs in the BF-HMW615 genome. The MAUVE Contig 
Mover module was used to order and orient the BF-HMW615 
contigs relative to the BF638R reference genome.'"' When needed, 
nucleotide or amino acid sequences were aligned using ClustalW. 

Annotation 

The sequences in GenBank were annotated by the Broad 
Institute as part of the genome sequencing but more than 80% 
of the proteins (122/144) were annotated as "hypothetical." 
To augment this annotation, we separately annotated all three 
sequences using the RAST Server, which assigns proteins to sub- 
systems. With the addition of the RAST information, 55 pro- 
teins remained annotated as hypothetical. Finally, we submitted 
the protein sequences to Phyre2-fold recognition analysis*^ and 
succeeded in obtaining a predicted protein function for all but 
17 proteins (38 were annotated at greater than 50% confidence 
and 25/55 at greater than 90% confidence) (Table S2). 

Mating of BF-HMW615 with B.fragilis 638R 

BF638R and BF-HMW615 were grown overnight anaero- 
bically at 37 °C in BHIS. 200 |xl of each culture were mixed, 
centrifuged at 13000 g for 2 min and then reconstituted in 100 
|jl1 of BHIS. The mixed cells were then spread on BHI plates 
(Anaerobe systems, CA, USA) and incubated anaerobically over- 
night at 37 °C. The growth on the plate was harvested and then 
resuspended (mating mix) in -1.5 ml 10% glycerol TSB storage 
media (Key Scientific Products Inc., TX, USA). A 150 |jl1 aliquot 
of mating mix was plated on BHIS plates containing tetracy- 
cline (1 (JLg/ml) and rifampicin (lOjJLg/ml). The frequency of the 
transfer was estimated to be 1 x 10 ' (number of transconjugants 
[BF638R-CtnHyb] /numbers of donors [BF HMW615]). 

Potential tetracycline-mediated increase in mating frequency 
was measured by three methods with slight alterations to the 



regular mating procedure. 1) BF-HMW615 was grown overnight 
in presence of tetracycline (1 |JLg/ml) prior to regular mating. 2) 
BF-HMW615 was treated with tetracycline (1 \y.%lm\) for 1 h at 
37 °C. The cells were then washed twice with equal volumes of 
BHIS and the regular mating procedure with BF638R was fol- 
lowed. 3) Overnight grown cells of BF638R and BF-HMW615 
mix (200 |jl1 of each) were each plated on BHIS plates contain- 
ing a very low concentration of tetracycline (0.001 |Jtg/ml), incu- 
bated overnight anaerobically at 37 °C, and then the regular 
mating procedure was followed. 

Mating of BF-HMW615 with E. coli 

200 |xl of overnight cuhures of BF-HMW615 and E. coli 
AGIOO cells were mixed and plated on BHI plates and incubated 
overnight at 37 °C under anaerobic conditions. The cells were 
then harvested and pooled in 1 ml of broth and 50 |xl of cells 
were plated on LB plate containing 40 |xg/ml kanamycin and 
incubated aerobically at 37 °C. 

Cloning of resistance genes contained within CTnHyb and 
introduction into E. coli 

HMPREF1204_02969 (a gene coding for the aminoglyco- 
side 3 '-phosphotransferase that confers kanamycin resistance) 
and HMPREF1204_02965 (a meJA-like efflux transporter gene) 
were PGR amplified from BF-HMW615, ligated into pSportI 
(Invitrogen) and introduced into E. coli DH5ol. Transformants 
were selected and purified on ampicillin containing LB plates 
and the presence of the introduced genes were confirmed by 
PGR and sequencing. 

MIC determinations 

MICs were determined using E-test technology (BioMerieux) 
according to the manufacturer's directions. An inoculum of 
one McFarland unit was used on a Brucella Blood Agar plate 
(Anaerobe Systems, Morgan Hill, CA). 

Transcription of genes contained within CTnHyb genes in 
BF-HMW615 

RNA was prepared from cells using the Qiagen RNAeasy 
kit (Qiagen, CA, USA) according to manufacturer's direc- 
tions. The total RNA was enriched for mRNA by removing 
the majority of rRNA using the Ambion Microbe Express Kit 
(Life Technologies). cDNA was prepared using the Invitrogen 
Superscript kit (Life Technologies). cDNA was quantified by 
RNA-Seq analysis (Otogenetics, Norcross, USA). RNA-Seq files 
were analyzed using the Lasergene Genomics Suite (DNASTAR, 
Inc., Madison, USA). 

Registration of CTnHyb 

CTnHyb was registered in the transposon registry (http:// 
www.ucl.ac.uk/eastman/research/departments/microbial-dis- 
eases/tn) as Transposon 6243. The complete sequence and anno- 
tation has been submitted to NCBI GenBank (Submission ID 
1725586, to be released June 10, 2014). 

Conclusion 

CTnHyb represents a mechanism for a single Bacteroides 
isolate to become a reservoir for a variety of resistance genes, 
both from Bacteroides and other species; these genes can then 
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be transferred both to other Bacteroides and unrelated bacteria. 
Identifying the factors that increase CTn accumulation within 
a strain as well as factors that increase transfer of CTns to other 
bacteria is critical information that could lead to therapeutic reg- 
imens against resistance dissemination. 
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